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ABSTRACT 

Supcrmassive black holes (BH) are powerful sources of energy that are already in place 
at very early epochs of the Universe (by z = 6). Using hydrodynamical simulations 
of the formation of a massive M v i r = 5 x 10 11 M halo by z = 6 (the most massive 
progenitor of a cluster of M V - 1T = 2 x 10 15 M Q at z = 0), we evaluate the impact 
of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and 
gas distribution inside this massive halo. We find that SN feedback has a marginal 
influence on the stellar structure, and no influence on the mass distribution on large 
scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge 
mass content by quenching star formation when the BH is self-regulating, and by 
depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH 
proceeds first by a rapid Eddington-limitcd period fed by direct cold filamentary infall. 
When the energy delivered by the AGN is sufficiently large to unbind the cold gas of 
the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow 
and clump migration through the galactic disc triggered by merger-induced torques. 
The feedback from the AGN has also a severe consequence on the baryon mass content 
within the halo, producing large-scale hot superwinds, able to blow away some of the 
cold filamentary material from the centre and reduce the baryon fraction by more 
than 30 per cent within the halo's virial radius. Thus in the very young universe, 
AGN feedback is likely to be a key process, shaping the properties of the most massive 
galaxies. 

Key words: cosmology: theory — galaxies: high redshift — galaxies: formation — 
galaxies: haloes — galaxies: active — methods: numerical 



1 INTRODUCTION 

Supermassive black holes (BH) are ubiquitous in the 
centres of galaxies with masses scaling with those of 
their hosts (Magorrian et al. 1998; Tremaine et al. 2002; 
Haring & Rix 2004), and with the most massive BHs found 
in the centres of the brightest cluster galaxies (BCGs). These 
observations suggest that a self-regulating process is at play 
in galaxy interiors involving some source of feedback en- 
ergy produced by accretion of material onto their central 
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BHs (Silk & Rees 1998; King 2003; Wyithe & Loeb 2003). 
The reality of such Active Galactic Nuclei (AGN) is con- 
firmed by direct detections of radio jets and X-ray cav- 
ities in the core of groups and clusters of galaxies (e.g. 
Boehringer et al. 1993), and also by observations of broad 
absorption lines in the spectra of quasars (e.g. Chartas et al. 
2003) 

Semi-analytical models of galaxy forma- 
tion (Croton et al. 2006; Bower et al. 2006; Cattaneo et al. 
2006; Somerville et al. 2008) with AGN feedback quench 
star formation in massive elliptical galaxies, and reproduce 
the bright end of the galaxy luminosity function with 
a sharp cut-off that is in good agreement with observa- 
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tions (Bell et al. 2003; Panter et al. 2007). Self-regulated 
growth of BHs with AGN feedback has also been modeled 
in various hydrodynamical cosmological simulations, repro- 
ducing the tight constraints on BH and galaxy relationships 
at z = 0, and confirming that this feedback strongly 
affects massive galaxy formation (Sijacki et al. 2007; 
Di Matteo et al. 2008; Booth & Schaye 2009; Dubois et al. 
2012; Kimm et al. 2012). These numerical models man- 
age to prevent the cooling catastrophe by regulating 
cooling flows in galaxy clusters, a necessary condition 
for obtaining massive galaxies with low star formation 
rates (Dubois et al. 2010, 2011; McCarthy et al. 2010, 2011; 
Teyssier et al. 2011). 

For scenarios without strong feedback, density profiles 
show large mass concentrations in the core regions of galaxy 
clusters that are at odds with observations. This is a conse- 
quence of the catastrophic cooling of the gas, which drives 
too many baryons into the cores of halos and consequently 
increases their gravitational potential wells via adiabatic 
contraction (Blumenthal et al. 1986; Gnedin et al. 2004), 
steepening their DM density profiles. This generic process 
has been suggested to occur early on in the mass assem- 
bly of clusters, and could be avoided by preventing baryons 
from cooling excessively and/or by blowing out large quan- 
tities of gas (Read & Gilmore 2005; Pontzen & Governato 
2012; Martizzi et al. 2012). As the gas in the intra-cluster 
medium is approximately in local hydrostatic equilibrium, 
gas thermodynamical properties are likely to be set up by 
the gas density distribution and the degree to which adi- 
abatic contraction has developed or has been prevented. 
Therefore without AGN feedback, gas in the intra-cluster 
medium is strongly depleted, giving unrealistic values for 
the temperature and entropy in cluster cores compared to 
X-ray observations (Tornatore et al. 2003; Valdarnini 2003; 
Nagai et al. 2007). The feedback activity at high redshift is 
likely to play an essential role in establishing and shaping 
some of the properties in clusters of galaxies, including the 
amount of gas present within the halo, the quantity of cold 
gas available to form stars, and the total stellar mass. 

Indeed, observations of distant and bright quasars sug- 
gest that strong feedback activity was present in the early 
epochs of galaxy formation before z = 6 (Fan et al. 2006; 
Jiang et al. 2009; Mortlock et al. 2011). Massive BHs with 
several 10 8-9 Mq are already in place, deeply embedded in 
dense gaseous regions, with rapid accretion onto them at 
rates close to the Eddington limit (Volonteri & Rees 2006; 
Willott et al. 2010; Treister et al. 2010) triggering powerful 
outflows. But, it is still very difficult to determine from high 
redshift observations the impact of early AGN feedback on 
the gas content of galaxies and the mass distribution within 
their halos. 

There are fundamental differences in the nature 
of the gas accreted onto galaxies at different epochs: 
red and massive galaxies are surrounded by hot and 
diffuse gas such as the intra-cluster medium de- 
tected in X-rays (Gursky et al. 1971; Forman et al. 1972; 
Markevitch et al. 1998; Vikhlinin et al. 2005; Pratt et al. 
2007), while theoretical and numerical models of ACDM 
universes predict that young and gas-rich galaxies 
obtain their gas through cold isothermal and col- 
limated flows (Binney 1977; Birnboim & Dekel 2003; 
Keres et al. 2005, 2009; Ocvirk et al. 2008; Brooks et al. 



2009; Dekel et al. 2009) that are still to be confirmed 
by direct detection (Kimm et al. 2011; Goerdt et al. 
2012). Though the interaction of AGN (radio jets) with 
their surrounding, hot intra-cluster gas has been ex- 
tensively studied (Omma et al. 2004; Ruszkowski et al. 
2004; Bruggen et al. 2005; Brighenti & Mathews 2006; 
Sternberg et al. 2007; Gaspari et al. 2011), little is known 
about how quasars at high redshift interact with cold in- 
falling streams of gas, and in particular whether these col- 
limated structures can survive the energy released through 
these powerful outflows. While feedback from supernova ex- 
plosions alone is not able to stop the streaming of cold 
flows towards galaxies (Powell et al. 2011), the larger en- 
ergies involved in AGN activity can potentially alter the 
motions and the morphology of the filamentary gas. In a 
large cosmological simulation box, Di Matteo et al. (2012) 
(see also Khandai et al. 2012) have shown that the primor- 
dial, most massive BHs can grow uninterrupted via cold fil- 
amentary infall, and despite AGN feedback, can acquire up 
to several billion solar masses at z = 6 in the most mas- 
sive halos (M v ir ~ 10 12_13 Mq) and inside the most active 
galaxies at this redshift, consistent with the observations of 
the most distant quasars (Wang et al. 2010, 2011). Similar 
conclusions have been reached previously by other authors 
using resimulations of individual massive objects (Li et al. 
2007; Sijacki et al. 2009). 

Using well-resolved hydrodynamical cosmological sim- 
ulations, Dubois et al. (2012) have shown that, in the ab- 
sence of feedback, the accretion onto the central bulge (and 
potentially a central BH) of massive galaxies is two-fold: 
during an early phase, cold streams feed directly the cen- 
tre of the galaxy; in a second phase, they contribute to the 
formation of proto-galactic discs by gaining angular momen- 
tum (Pichon et al. 2011). A large amount of gas is accreted 
onto the disc, maintains it to be gravitationally unstable 
and will rapidly feed the compact bulge through migration 
of clumps (see also Bournaud et al. 2011). 

We will address here the question of massive galaxy for- 
mation at high redshift from another angle while focusing on 
the impact of AGN feedback onto the early assembly phase 
of a progenitor of a massive cluster of galaxies - but not nec- 
essarily the most massive object at z = 6 - using sufficient 
resolution to capture the detailed structure of the galaxy and 
of the inflow/outflow onto it (Ax — 15 pc). We will ques- 
tion whether feedback from AGN is able to rapidly quench 
star formation inside the central galaxy, and what kind of 
outflows can form from the central source of energy. Our 
model of AGN feedback has been calibrated to reproduce 
observational constraints at z — (such as BH mass den- 
sity and Msn-Mh, Dubois et al. 2012) using kpc-resolution 
hydrodynamical cosmological simulations, and no recalibra- 
tion of the model is performed in what follows, unless stated 
explicitly. Standard recipes are used to describe the physics 
of star formation and stellar feedback. Hence in this paper 
we investigate the consequences of such a "canonical" feed- 
back model on the early galaxy evolution and its associated 
inflow/outflow. 

The paper is organised as follows: in section 2, we 
present the details of the numerical implementation. We 
report the impact of AGN feedback on the galaxy prop- 
erties (section 3), the accretion rate onto the central BH 
(section 4), and the impact on the gas content within the 



AGN and high-redshift galaxy formation 3 



halo (section 5). Finally we discuss the implications of such 
strong AGN feedback (section 6). 



2 NUMERICAL SET-UP 

The initial conditions are designed to yield a rare and mas- 
sive cluster halo with a virial mass of M v i r (z = 0) = 
2 x 10 15 M at z = 0. We assume a ACDM cosmol- 
ogy with total matter density Q, m = 0.27, baryon density 
Qb = 0.045, dark energy density S1a = 0.73, amplitude 
of the matter power spectrum erg = 0.8 and Hubble con- 
stant Ho = 70 km s" 1 Mpc" 1 consistent with the WMAP 
7-year data (Komatsu et al. 2011). The box size of our sim- 
ulations is Lbox = 100 ft" 1 Mpc, within which we define a 
high resolution region where high resolution DM particles 
(M res = 9 x 10 5 ft" 1 M Q ) fill a volume of size 2 r v i r around 
the main progenitor of the targeted massive cluster halo at 
z = 6. The mass of DM particles is slowly downgraded from 
inside-out. The mass high resolution region fills a sphere of 
2.3 ft" 1 Mpc comoving radius in the initial conditions, then 
lower resolution DM particles fill concentric shells as fol- 
low: 7 x 10 6 ft" 1 M from 2.3-2.8 ft" 1 Mpc, 6 x 10 7 ft" 1 M Q 
from 2.8-3.5 ft" 1 Mpc, 5 x 10 s ft" 1 M from 3.5-5 ft" 1 Mpc, 
4 x 10 9 ft" 1 M from 5-7.5 ft" 1 Mpc, 3 x 10 10 ft" 1 M Q from 
7.5-15 ft" 1 Mpc, and 2 x 10 11 ft" 1 M filling the rest of the 
box. To smooth out the contribution of DM particles to the 
gravitational potential in highly refined regions of the gas 
flows, their mass is projected up to the level 16 of the grid 
(while the grid reaches up level 20 at z = 6) for the compu- 
tation of the gravitational potential that correspond to an 
effective softening length of 300 pc at redshift 6. This is 20 
times the softening length for the hydro solver at this red- 
shift. The progenitor of the halo has a mass and virial radius 
at z = 6 of M v ir(z = 6) = 5 x 10 11 M and r vir = 35 kpc. 

The simulations are run with the Adaptive Mesh Refine- 
ment code RAMSES (Teyssier 2002). The evolution of the gas 
is followed using a second-order unsplit Godunov scheme for 
the Euler equations. The HLLC Riemann solver with a first- 
order MinMod Total Variation Diminishing scheme to re- 
construct the interpolated variables from their cell-centered 
values is used to compute fluxes at cell interfaces. Collision- 
less particles (DM and star particles) are evolved using a 
particle-mesh solver with a Cloud-In-Cell interpolation. The 
mesh is refined up to Aa; = 15 pc using a quasi-Lagrangian 
strategy: when more than 8 DM particles lie in a cell, or 
if the baryon density is larger than 8 times the initial DM 
resolution. 

Gas is allowed to cool by H and He cooling with 
a contribution from metals using a Sutherland & Dopita 
(1993) model for temperatures above 10 4 K, and metal fine- 
structure cooling below 10 4 K following Rosen & Bregman 

(1995) . The gas temperature is not allowed to fall below 
To = 100 K. Heating from a uniform UV background takes 
place after redshift z re ion = 8.5 following Haardt & Madau 

(1996) . Metallicity is modeled as a passive variable for the 
gas (whose composition is assumed to be solar) and is al- 
tered by the injection of gas ejecta during supernovae (SNe) 
explosions, assuming a constant yield of 0.1. We assume an 
initial metallicity of Z = 10~ 3 Z . 

The star formation process is modeled with a Schmidt 
law: p* = e*p/tff , where p* is the star formation rate den- 



sity, e* the constant star formation efficiency, and ta the 
local free-fall time of the gas. We choose a low star for- 
mation efficiency e» = 0.01 consistent with observations of 
giant molecular clouds (Krumholz & Tan 2007). Star for- 
mation is allowed in regions exceeding a gas density thresh- 
old of po = 50 Hem" 3 . The gas pressure is artificially en- 
hanced above p > po assuming a polytropic equation of state 
T = To(p/po) re_1 with polytropic index k = 2 to avoid exces- 
sive gas fragmentation. Feedback from SNe releases 10 51 erg 
per 10 M assuming a Salpeter IMF where 77SN = 0.1 frac- 
tion of the IMF is composed of type II core-collapse SNe. 
The energy from SNe is coupled to the gas with a kinetic im- 
plementation from Dubois & Teyssier (2008), where mass, 
momentum and kinetic energy profiles are imposed to mimic 
a Sedov blast wave. Neither feedback from stellar winds nor 
type la SNe are taken into account. 

We use the same "canonical" AGN feedback modeling 
employed in Dubois et al. (2012). BHs are created at loci 
where the gas density is larger than the density thresh- 
old for star formation po with an initial seed mass of 
10 5 M . We assume that very massive BHs form from the 
direct collapse of the low-angular momentum gas in pri- 
mordial halos (Begelman et al. 2006). In order to avoid the 
formation of multiple BHs in the same galaxy, BHs are 
not allowed to form at distances smaller than 5 kpc from 
any other BH particle. The accretion rate onto BHs fol- 
lows the Bondi-Hoyle-Lyttleton (Bondi 1952) rate Mbh = 
4naG 2 Mi ll p/(c 2 s + u 2 ) 3/2 , where M B h is the BH mass, p 
is the average gas density, c s is the average sound speed, u 
is the average gas velocity relative to the BH velocity, and 
a is a dimensionless boost factor with a = (p/po) 2 when 
p > po and a = 1 otherwise (Booth & Schaye 2009) in or- 
der to account for our inability to capture the colder and 
higher density regions of the ISM. The effective accretion 
rate onto BHs is capped at the Eddington accretion rate: 
A^Edd = 47rGJ\fBH?rip/(erO"Tc), where ot is the Thompson 
cross-section, c is the speed of light, m p is the proton mass, 
and e r is the radiative efficiency, assumed to be equal to 
e r = 0.1 for the Shakura & Sunyaev (1973) accretion onto 
a Schwarzschild BH. In order to avoid spurious oscillations 
of the BH in the gravitational potential well due to exter- 
nal perturbations and finite resolution effects, we introduce 
a drag force that mimics the dynamical friction exerted by 
the gas onto a massive particle. This dynamical friction is 
proportional to F D f = f zas A-Kap(GMBn/c s ) 2 , where / gas is 
a fudge factor whose value is between and 2 and is a func- 
tion of the mach number M = u/c, < 1 (Ostriker 1999; 
Chapon et al. 2011), and where we introduce the boost fac- 
tor a for the same reasons than stated above. 

The AGN feedback is a combination of two differ- 
ent modes, the so-called radio mode operating when \ = 
AfBH/MEdd < 0.01 and the quasar mode active otherwise. 
The quasar mode corresponds to an isotropic injection of 
thermal energy into the gas within a sphere of radius Aa;, 
at an energy deposition rate: -Eagn = e{e T MBnc , where 
6f = 0.15 for the quasar mode is a free parameter cho- 
sen to reproduce the Mbh-My,, MeH-Cb, and BH density 
in our local Universe (see Dubois et al. 2012). At low accre- 
tion rates on the other hand, the radio mode deposits the 
AGN feedback energy into a bipolar outflow with a jet veloc- 
ity of 10 4 km/s into a cylinder with a cross-section of radius 
Aa; and height 2 Aa; following Omma et al. (2004) (more de- 
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Figure 1. Images of the gas density (top), and the stellar luminosity in the SDSS ugr filter bands (bottom) of the central galaxy seen 
face on at z = 6 for the nofeedback case (left), the noAGN case (middle), and the AGN case (right). The image sizes are 5 kpc length 
and deep. The nofeedback / no AGN cases have comparable gas densities, stellar intensities and number of clumps, while the AGN case 
is different. 



tails about the jet implementation are given in Dubois et al. 
2010). The efficiency of the radio mode is larger with et — 1. 
We insist that we are taking the parameters of the AGN 
model at face- value (as we do for the parameters of star for- 
mation and SN feedback) calibrated on lower resolution hy- 
dro simulations that reproduce the Magorrian et al. (1998) 
relations at z = 0. 

We are using three sets of simulations to asses the role 
of the different modes of feedback, focusing on the role of 
AGN feedback on galactic evolution and the mass content 
in the halo. 

• nofeedback: the reference simulation does not include 
any type of feedback. Only gas cooling, background UV 
heating and star formation are allowed. 

• no AGN: the second simulation includes all the physics 
of the nofeedback run in addition to feedback from type II 
SNe, but no AGN feedback. 

• AGN: the final simulation includes all the physics of 
the no AGN run and feedback from AGN in a radio/quasar 
mode as described above. 

We also complemented these runs with two other sim- 
ulations which include AGN feedback with different flavors: 
one where only the radio (jet) mode is active to test the 
effect of the anisotropy of the energy injection, and another 
with the dual radio/quasar mode but with the AGN effi- 
ciency ef reduced by a factor of 15 to test how this affects 
the energy released by BHs. We do not detail the results from 
these simulation runs, but will punctually use them to rein- 



force the robustness of our findings concerning the impact 
of AGN feedback on structure formation (see Appendices A 
and B). 

We performed some additional simulations to test the 
impact of changing the resolution in DM mass and in min- 
imum cell size on our final results. We find no significant 
difference with the results obtained using our canonical res- 
olution, even though the resolution elements are changed by 
one order of magnitude (in DM mass and minimum cell size, 
see Appendix C for further details) . 



3 QUENCHING STAR FORMATION AND 
CATASTROPHIC GAS COOLING 

Fig. 1 shows the gas density and the stellar luminosity 
(SDSS ugr filters) images of the central galaxy at z = 6 
for the three different runs, revealing in each case, a massive 
bright bulge component with multiple clumps in rotation 
with the disc component. However the nofeedback and the 
noAGN simulations do not show large differences: they share 
a similar number of gas/stellar clumps with similar densi- 
ties/intensities. In contrast, as a consequence of the impact 
of AGN feedback on the baryon content of the galaxy, the 
AGN case has fewer clumps and less dense gas available in 
the galaxy. We also note that recent AGN activity in the 
bulge creates a hole in the gas distribution in the centre of 
the galaxy in the AGN run. 

The measured star formation rates (SFRs) in the galaxy 
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Figure 2. Top: SFR history using the stars contained within 
0.1r v ; r for the nofeedback case (black), the no AGN case (blue), 
and for the AGN case (red). Bottom: SFR within r^, = 50 pc 
(same color coding as top panel), and BH accretion rate (grey 
solid), Eddington accretion rate (grey dashed), and 1 per cent 
of Eddington accretion rate (grey dotted). Merger events for 
the AGN run are marked with vertical red lines. SFRs of the 
nofeedback /no AGN cases show almost no difference, while the 
SFR in the AGN case is significantly reduced when the central 
BH exhibits a departure from its Eddington-limited growth after 




redshift 

Figure 3. Mass accretion rate at 0.1r v ; r as a function of 
redshift for the nofeedback case (black), no AGN case (blue), 
AGN case (red). Smooth accretion onto the galaxy is reduced 
in the AGN case after z = 8, but mergers continue to proceed. 



for the stars enclosed within 0.1r v i r at z = 6, indicate that, 
without AGN, the galaxy maintains high levels of star for- 
mation (> lOOMoyr -1 ) (see Fig. 2). Feedback from SNe 
has only a marginal effect on the amount of stars formed in 
the galaxy. This is also the case if the SFR measurement is 
done for the stars in the bulge of the galaxy (r < where 
i"h = 50 pc). The bulge and disc components are decomposed 
by fitting two exponentially decreasing profiles £i(r) = 
£i,o exp(— r/n) 1 plus one gaussian profile (for the bar com- 
ponent if any) on the stellar surface densities. Masses of the 
bulge and of the disc are defined as Mi = Ei^^rf . The mea- 
sured bulge stellar mass is the same for the nofeedback run 
and the no AGN run (Mb.nof.noA = 1.9 X 10 10 M Q at z = 6), 
while the disc stellar mass is reduced by 20 per cent due 
to the presence of SNe (from M d , nof = 3.0 x 10 10 M Q to 
Md,noA = 2.4 x 1O 1O M ). 

As already suggested qualitatively by Fig. 1, the SFR 
measurement for the AGN run confirms that AGN feedback 
has a dramatic impact on the amount of stars formed in 
the galaxy. The global SFR is reduced to a few 10 Mq yr~\ 
and the SFR in the bulge is almost completely quenched 
to less than 1 Mq yr~ , two orders of magnitude below the 
no AGN case. Indeed, it has severe consequences for the 
bulge stellar mass that is reduced by a factor 3 with Mb, a = 
6.2 x 10 9 M at z = 6. The bulge, despite the impact of AGN 
feedback, is more compact than the most compact bulges of 
comparable stellar mass at 1 < z < 3 (Bruce et al. 2012), 
indicating that bulges follow a strong evolution with redshift 
and are more compact in the distant Universe (Daddi et al. 
2005; Khochfar & Silk 2006; Cimatti et al. 2008). The disc 
stellar mass is also reduced to M^a = 1-2 x 1O 1O M0, 
mainly because the AGN clears the centre of the halo, caus- 
ing the accretion of new material to proceed more quietly, 
while for the no AGN run, energy from SNe is not suffi- 
cient to significantly impact the accretion rate at 0.1r v i r 
(see Fig. 3). Indeed after 2 = 8 the mass inflow onto 
the galaxy proceeds with a smooth component, correspond- 
ing to the accretion of cold streams, at a few lOMQyr -1 
for the AGN run, as opposed to > lOOMQyr -1 for the 
nofeedback I no AGN simulation. Bursty episodes mark the 
accretion of small satellites. In summary, AGN are able to 
substantially reduce the smooth part of the accretion flow, 
but the bursty episodes of accretion persist. 

The consequences of AGN feedback on the mass pro- 
files of baryons in the galaxy are also dramatic (Fig. 4). 
The stellar and gas densities are reduced both in the bulge 
and the disc with the formation of a gap in the gas dis- 
tribution within the bulge due to the outflow produced by 
the central BH. Again, the strongest impact of the AGN is 
on the bulge component, but the disc is also significantly 
altered because the accretion onto the galaxy is quenched 
efficiently. A natural outcome is that the peak of the cir- 
cular velocity is reduced from 900 to 500 kms - , together 
with a reduced velocity dispersion of the stars within the 
bulge from 450 to 250 kms -1 , consistent with velocity dis- 



1 Bulges are more usually defined by a Sersic profile with an index 
depending on its morphological classification (bulge or pseudo- 
bulge). We do the simple assumption that our simulated stellar 
bulge profiles can be described by a Sersic profile and index of 1, 
corresponding to an exponentially decreasing profile. 




Figure 4. Density profiles of all the components (gas+stars+DM) (solid lines), stars (left panel dashed), gas (middle panel dashed), 
and DM (right panel dashed) for the nofeedback case (black), the noAGN case (blue), and the AGN case (red) at z = 6. The vertical 
grey solid line indicates the resolution limit, the vertical grey dashed lines indicates the scale radius of the bulge, and the vertical grey 
dotted line the scale radius of the disc. 



persions observed in z — 6 galaxies inferred from velocity 
widths of CO lines (Wang et al. 2010). Note also the deficit 
of gas between the disc radius at r — 2 kpc and 10 kpc that 
is due to efficient ejection of gas around the galaxy. The fact 
that the stellar and gas distribution in the no AGN case 
are similar to their nofeedback counterparts, even in the 
disc component, means that SNe do not efficiently remove 
gas or quench the gas infall onto the galaxy, and that with 
similar gas content in the galaxy, identical stellar distribu- 
tions are found. The ability of AGN feedback to prevent 
the baryons from collapsing and/or to remove the gas from 
the galaxy early-on weakens the adiabatic contraction of 
the dark matter particles (Blumenthal et al. 1986) and is 
therefore key to understanding mass density distributions 
in galaxy clusters (Dubois et al. 2010, 2011; Teyssier et al. 
2011; Martizzi et al. 2012). Indeed, we see that with AGN 
feedback, the density profiles of gas, stars and DM are shal- 
lower inside the galaxy, as AGN feedback prevents the cold 
gas from collapsing too much. The measured mass of DM in 
the centre of the halo within 0.1 r v i r has been reduced by 35 
percenttoM D M,A(0.1r vh .) = 3.8xl0 10 M Q in the A GN case 
compared to the mass M DM ,„of(0.1r vh .) = 6.0 x 1O 1O M 
in the nofeedback case at z = 6 (Mdm,iioa(0.1 r V ir) — 
5.4 x 10 10 M in the no AGN case). 



4 WHAT DRIVES THE AGN ACTIVITY? 

According to Fig. 2, up to 600 Myrs after the Big Bang 
the SFRs for the AGN case and the no AGN case are com- 
parable. The onset of SFR quenching at z ~ 7.5 coincides 
with the redshift when the BH stops accreting gas at the 
Eddington limit (see bottom panel of Fig. 2). After z = 8, 
the growth of the central BH is regulated by its own quasar- 
mode feedback which can efficiently remove gas from the 
galaxy. The strong outflows produced by the central quasar 
prevent gas accretion onto the galaxy (see Fig. 3), thereby 
starving the gas reservoir available for star formation. The 
accretion of gas for the AGN case looks spikier compared 
to the no AGN case or nofeedback case due to efficient re- 
moval of the diffuse component surrounding the galaxy with 
low ram-pressure, while the dense star forming gas contin- 
ues to accrete onto the galaxy (see also Fig. 10). Because 



BH accretion rates hover at around a tenth of the Edding- 
ton rate, the AGN feedback proceeds in quasar mode all 
the way down to z = 6. Recall that in this mode energy 
is deposited isotropically. However due to the presence of a 
dense disc component surrounding the bulge of the galaxy, 
the AGN launches bipolar outflows with large opening an- 
gles from the galaxy. 

To test the importance of the directionality and colli- 
mation of the outflow on our results, we run a similar sim- 
ulation to the AGN case where only the radio (jet) mode 
is allowed. Despite the absence of opening angle in the jet 
input, and its small size input (Ax = 15 pc) compared to 
the extent of the galactic disc (2 kpc), we reach similar con- 
clusions because the jet propagates as a light jet into the 
hot diffuse gas, and almost immediately spreads into a wide 
outflow (see Appendix A). 

Analysing the halo in the nofeedback simula- 
tion, Dubois et al. (2012) have shown that the cold 
flows feed the central bulge directly until z ~ 8. After this 
redshift, the filaments wrap around the galaxy because they 
gain angular momentum from the tidal field of the large 
scale cosmic structures (Pichon et al. 2011). This gives rise 
to the formation of a large disc component surrounding 
the bulge that is gravitationally unstable (with Toomre 
parameter of Q ~ 1) causing disc instabilities and clump 
migration (Elmegreen et al. 2008; Bournaud et al. 2011) to 
bring large quantities of gas towards the bulge. Interestingly, 
the central BH stops accreting at the Eddington rate at the 
moment when the gas accretion onto the bulge transitions 
from direct infall from cold streams to the feeding through 
disc instabilities (see Fig. 10 in Dubois et al. 2012). This 
suggests that the cold flows are directly responsible for 
maintaining sufficiently large amounts of gas in the bulge, 
so as to allow the BH to accrete at its maximum rate. 

Several mergers that correlate with short bursts of star 
formation can be identified in the SFR evolution (vertical 
red lines in Fig. 2). The response of the SFR to such mergers 
is different depending on the presence or absence of AGN. In 
the nofeedback and noAGN cases, the SFR after a merging 
event reverts to the SFR levels present before the merger. 
Thus, the outflow produced by SNe during the starburst is 
not strong enough to regulate the gas content in the galaxy. 
In contrast, as soon as the BH enters a sub- Eddington mode 
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Figure 5. Time sequence (as labelled) around a peak of starburst and quasar activity triggered by a merger of 1:4 ratio at z ~ 6.5 of 
the gas density of the central galaxy. Red circles correspond to the position of the two most massive BHs. The projections of the images 
are 9 kpc deep. The panels on top of each image show the bolometric luminosity of the central BH in Lq (black) with the SFR (green) 
of the central galaxy and the scaled mass accretion rate onto the bulge measured at 80 pc (blue). The SFR and mass accretion rate 
values in Mq yr — 1 are divided by a factor 10 and 30 respectively for a better comparison to the luminosity of the central BH. The burst 
of SFR happens at the passage of the satellite at the apocentre. The rapid quasar activity is triggered by the enhanced migration of a 
massive clump of dense gas (black circles) whose motion is perturbed by the presence of the satellite. Note that the peak of SFR and 
AGN activity are not synchronized. The final merger between the two galaxies is complete at z = 6.2 . 



of accretion (after z ~ 8) in the AGN case, mergers provoke 
a large drop of the SFR in the bulge, following the star- 
burst: the SFR decreases by one order of magnitude below 
its pre-merger value. This highlights the fact that mergers 
trigger large inflows of gas towards the centre of the merging 
galaxies and can efficiently feed the BHs with fresh gas that 
activates a strong period of AGN activity able to blow away 
the surrounding gas (Di Matteo et al. 2005; Springel et al. 
2005; Hopkins et al. 2006). 

The largest galaxy merger in the simulation occurring 
at z ~ 6.5 with mass ratio 1:4 is illustrated in Fig. 5 with 
a time sequence of the gas density along with the central 
quasar bolometric luminosity and SFR of the central galaxy. 
We see that the peak of the quasar activity at z = 6.52 is 
not synchronized with the peak of the starburst at z = 6.61, 



corresponding to a time delay of 15 Myrs equal to the dy- 
namical time of the galaxy td — 17 Myr (computed for the 
enclosed mass within 0.1r v i r ). The latter occurs when the 
satellite passes the apocentre, while the former is correlated 
with the sharp increase of the mass inflow rate to the bulge 
when a satellite exerts a global torque on the central galaxy. 
An important point is that it is not because the two galaxies 
merge that the AGN activity peaks (the merger of the two 
galaxies is not complete by z = 6.2), but rather because the 
satellite applies this strong torque on the disc of the cen- 
tral galaxy. In response, the gas looses angular momentum 
and plunges radially towards the bulge, and clumps of gas 
undergo a faster migration towards the bulge. This effect is 
illustrated with the clump circled in black in Fig. 5 with gas 
mass M gi c = 5 x 10 8 Mq. The final capture of the clump 
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Figure 6. Top: BH mass as a function of redshift. Bottom: Feed- 
back luminosity from type II SNe for the noAGN case (blue) and 
for the AGN case (red), and the feedback luminosity from the 
central AGN in the AGN run (grey). The energy output from 
the central AGN largely dominates the total energy budget from 
feedback. 



corresponds to the burst in the bulge accretion rate followed 
by the peak of the AGN activity. The idealized picture de- 
scribed in Bournaud et al. (2011) where the accretion onto 
the central BH is triggered by the rapid migration of gas 
clumps into the bulge should therefore be supplemented by 
the following: external perturbations exert strong torques on 
the disc of the galaxy that further accelerate the accretion 
of clumps onto bulges. 

Computing the energy budget involved in the differ- 
ent feedback processes helps us to explain why AGN feed- 
back is able to unbind the gas from the galactic cen- 
tre while SNe are extremely inefficient. The total energy 
injected by the central BH through the quasar mode is 
-Eagn = e/e r MeHC 2 ~ 2 x 10 60 erg where we take Mbh to 
be the mass of the central BH Mbh = 8x 10 7 M at z — 6 
(Fig. 6), whereas the energy liberated by SN explosions oc- 
curring in the bulge is -Esn.a = 77sNM b , A (lO 51 erg/lOM ) ~ 
6 x 10 58 erg for the AGN case (£ S N, m A ~ 2 x 10 59 erg 
for the no AGN case). For comparison, the binding energy 
of the bulge is £ b , A ^ GMljr h ~ 6 x 10 58 erg for the 
AGN run (£b,noA ~ 6 x 10 59 erg for the noAGN run). 
Hence the energy from SNe, in the noAGN case is com- 
parable to the binding energy of the bulge but because a 
non-negligible fraction of the SNe energy is lost through ra- 
diative cooling, SNe cannot remove the gas from the deep 
gravitational potential well. AGN feedback on the other 
hand delivers an amount of energy larger by more than one 
order of magnitude than the total binding energy of the 
bulge, and cleans out its gas content. Fig. 6 shows that the 
AGN feedback luminosity Lagn = -Eagn after z — 9 is al- 
ways larger than the feedback luminosity from type II SNe 



Lsnii = SFR x 77SN(10 51 erg/10M Q ). Also the luminosity of 
the AGN needs to reach high values (Lagn = 10 Lq at 
2 = 7.6) before it starts to self-regulate the BH growth: this 
is the point at which the AGN finally unbinds the gas within 
the bulge. 

The final BH mass at z = 6 is Mbh = 8 x 1O 7 M , 
somewhat low compared to the brightest quasars observed 
at z = 6 (Fan et al. 2006). Although it is not the aim 
of this paper to match these observations (we are inter- 
ested in the early evolution of the progenitor of a typi- 
cal supermassive cluster), we note that the mass of the 
halo probed in that simulation turned out to be probably 
below the mass needed to host these exceptional BHs if 
"canonical" AGN feedback is taken into account. First, the 
typical SFRs observed in the far-infrared suggest that the 
brightest quasars produce stars at lO 3 M yr _1 , somewhat 
higher than even the star formation in our nofeedback sim- 
ulation (SFR~ 100 - 2OOM yr _1 ). Secondly, the maxi- 
mum halo masses obtained for standard ACDM and for 
which the SFR is potentially the highest are of the or- 
der of several 10 12 M 2 . The number density of observed 
quasars at z = 6 is extremely low with one detectable ob- 
ject per comoving Gpc 3 (Fan et al. 2001), which suggests 
that only the most extreme and rarest environments can 
be the hosts of the brightest quasars in the Universe. If we 
depart from the canonical model and choose a lower AGN 
efficiency, the BH has to grow to a larger mass in order to 
release the same amount of energy to reach self-regulation 
(as already demonstrated in Dubois et al. (2012) using a set 
of large cosmological simulation boxes with multiple halo 
masses). Our test run with an AGN feedback efficiency tt 
reduced by a factor of 15 results in a more massive BH of 
mass Mbh = 5 x 10 s M at z = 6.7. Note that this re- 
duced efficiency does not significantly change the impact of 
the AGN feedback on the SFR of the central galaxy, be- 
cause similar energy is released by the central BH (see Ap- 
pendix B). The BH mass is a factor four below the local 
MeH-ffb relationship (Tremaine et al. 2002), and a factor 
twelve above the local Mbh-M^ relationship (Haring & Rix 
2004) for our canonical model of AGN feedback. Observa- 
tions from distant galaxies suggest that such an evolving 
relationship MBH-Mb (positive evolution with redshift) is 
present (Decarli et al. 2010; Merloni et al. 2010), also con- 
firmed by cosmological simulations (Di Matteo et al. 2008; 
Booth & Schaye 2011; Dubois et al. 2012). 



5 REMOVING GAS FROM THE HALO 

AGN feedback has an important impact on the mass con- 
tent within the galaxy, and also on larger scales by virtue 
of it removing gas efficiently from the halo. Fig. 7 shows 
the gas density and gas temperature distribution in the 
halo at z — 6 for the three different runs. From a compar- 
ison of the nofeedback (upper panels) and no AGN (mid- 
dle panels) cases, it appears that SN explosions in the 

2 The maximum halo mass at z = 6 in the MassiveBH simulation 
from Di Matteo et al. (2012) is Af vir = 5 X 10 12 M . In addition, 
they find a BH mass of Mbh — 2 x 10 7 M for a typical average 
halo mass of M vir = 5 x 10 11 M (see fi gure 1 of DeGraf et al. 
2012) which is the halo mass simulated here. 
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Figure 7. Mass- weighted projected number densities (left panels), and temperatures (right panels), for the nofeedback case (top panels), 
the noAGN case (middle panels), and the AGN case (bottom panels) at z = 6. The projections are 80 kpc deep. The small and large 
circles are respectively 0.1r v j r and r v i r . Both simulations without AGN feedback show strikingly similar features with cold streams of gas 
connecting to the central galaxy in a continuous manner. In the AGN case, these streams are strongly disturbed and the gas temperature 
is increased and expands to larger distance. 
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Figure 8. Baryon fraction (solid), stellar fractions (dashed), 
gas fractions (dotted) at z = 6 for the nofeedback case (black), 
noAGN case (blue), and for the AGN case (red). The horizon- 
tal grey line is the universal baryon fraction. Baryon fractions 
in the nofeedback / no AGN cases at r v ; r (vertical grey line) are at 
the universal baryon fraction. The AGN feedback significantly re- 
duces the baryon fraction at r v ; r by 30 per cent below its universal 
value. 



no AGN run do not perturb the cold gas streaming toward 
the central galaxy: the same features are present as in the 
nofeedback case, and the cold gas connects all the way down 
to the central galaxy. The most obvious difference between 
the two simulations is that the hot gas filling the halo is 
slightly more extended in the simulation with SNe explo- 
sions. The difference with the AGN case however is strik- 
ing: very hot (T~ 10 7 K) gas pours into the circum-galactic 
medium and the dense cold streams are morphologically per- 
turbed within the halo. We will now quantify the extent to 
which AGN feedback is able to remove gas from the halo. 

Fig. 8 displays the baryon fraction as a function of ra- 
dius for the three simulations at z = 6. While we see that 
SN feedback alone has almost no effect on the total baryon 
fraction at the virial radius r v i r , the addition of AGN feed- 
back drives more than 30 per cent of the baryons out of 
the halo (beyond two 2 r v j r that is the size of the high res- 
olution region) and decreases the total stellar fraction (or 
mass) by a factor 2. This confirms that strong and early 
feedback activity can significantly reduce the total amount 
of baryons in the halo. Feedback at later times, might be 
less efficient because halos tend to be more spherical, mak- 
ing it more difficult for gas to escape from collapsed struc- 
tures (Peirani et al. 2012). We also demonstrate that a stan- 
dard recipe for SN feedback does not suffice to alter the 
baryonic mass content in a progenitor of a massive halo, 
whereas AGNs, with their large energy release (more than 
a factor 10 compared to the energy released by SNe), are 
good candidates for powerful early pre-heating. 

It is important to investigate the ability of AGN feed- 
back to efficiently remove the gas from the halo, and whether 
the different gas phases are able to survive within the hot 
atmosphere produced by ejecta from the central BH. We de- 
vide the gas into four different phases: the star forming gas 
(ISM) with gas density above p ^ pa, the collimated cold 
flows (F) with gas density 100p a v g ^ p < po and tempera- 
ture T < 10 5 K, the cold diffuse gas (CD) with p < 100p avg 
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Figure 9. Temperature-density diagrams of all the gas within 
r v ; r for the nofeedback case (top), no AGN case (middle) and the 
AGN case (bottom) at z = 6. Colors represent the logarithm 
of the metallicity of the gas while the black solid contours are 
the mass-weighted contours with 1 dex difference between each 
contour. The solid horizontal and vertical lines separate the four 
different phases of the gas (see text for details). The dashed hor- 
izontal line is the virial temperature of the halo at z = 6. 
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and T < 10 5 K, and the hot diffuse gas (HD) with p < po 
and T ^ 10 s K, where we define p avc as the cosmic average 
gas density. Fig. 9 shows temperature-density diagrams of 
the gas metallicity and of the gas mass distribution of the 
gas within the virial radius of the halo at z — 6. Even with- 
out SNe or AGN feedback, a hot phase with a temperature 
close to the virial temperature of the halo, T V - 1T = 6 x 10 6 K 
is present. This hot phase originates from cosmic gas that 
shock-heats as it accretes onto the halo, and for which the 
radiative cooling is not efficient enough (Birnboim & Dekel 
2003). When SNe deliver thermal energy to the gas, the tem- 
perature of the hot gas component rises above the virial tem- 
perature of the halo. AGN feedback increases the tempera- 
ture of the hot gas even further: the mass- weighted tempera- 
ture of the hot gas phase is now Thd,agn = 8 x 10 6 K instead 
of Thd.iioAGn = 3.5 x 10 6 K (similar to the nofeedback case). 
The bursts of AGN feedback also produce more mixing of 
metals in the hot diffuse g well as in the cold dif- 

fuse gas: we will confirm this effect below by analyzing the 
morphological structure of the filaments. Pushing the metal 
rich gas out of the halo centre and increasing its mixing 
are aspects of AGN feedback believed to play important 
roles in getting realistic flat metallicity profiles at late times 
in galaxy clusters (Leccardi & Molendi 2008; Fabjan et al. 
2010; McCarthy et al. 2010) 

Fig. 10 shows the mass fluxes measured at r v i r for the 
no AGN and AGN runs as a function of redshift. Outflow 
and inflow components (positive and negative fluxes) are 
separated according to their radial velocities: positive ra- 
dial velocities correspond to outflows and vice versa. The 
no AGN case reveals that the majority of the mass is ac- 
creted in the form of cold flows, with an increasing amount 
of inflowing diffuse (hot and cold) gas as the halo be- 
comes sufficiently massive to enter a phase where both 
cold flows and shock- heated gas coexist (Ocvirk et al. 2008). 
In the AGN run, cold diffuse accretion is quenched, hot 
diffuse accretion is significantly reduced, cold accretion is 
slightly reduced, and large amounts of hot gas are driven 
beyond the virial radius (> lOOMQyr" 1 ). The hot out- 
flowing phase manages to reduce the net flux to less than 
< lOOMQyr -1 entering the halo at z — 6 as opposed to 
a flux of ~ 250M Q yr~ 1 without AGN. The intensity of 
the mass outflow generated by the central AGN is com- 
parable to the SFR measured in the absence of feedback 
(~ lOOMQyr -1 , see Fig. 2), thus, the mass loading factor 
of the wind is large enough to self-regulate the SFR inside 
the galaxy. The outflow velocity of the hot gas at r v ; r is 
~ 500 km s -1 averaged over time during the removal of the 
hot gas from the halo, with values close to the galaxy that 
can reach several 1000 km s^ 1 consistent with outflow veloc- 
ities observed in powerful distant quasars (Maiolino et al. 
2005, 2012; Cano-Di'az et al. 2012). This value is larger than 
the escape velocity of the halo v csc — 350 km s -1 at z — 6 
and explains why the baryon fraction within the halo is effi- 
ciently reduced. Note that the outflow rate of the hot diffuse 
phase does not continuously increase but instead endures a 
sharp drop after z — 6.5 corresponding to the global de- 
cline of the central quasar activity after it enters its sub- 
Eddington phase (at z — 7.5, see Fig. 6) with a time lag. At 
late times, around z = 6, some of the cold filamentary gas 
is entrained by the hot gas and flows out beyond r v - a , but 
most of the gas continues to be accreted. 



Fig. 11 shows the integrated mass of star forming gas, 
cold filaments and hot gas at z = 6 for the three different 
simulations. SN feedback is not able to remove the cold star 
forming ISM from the central galaxy, and the filamentary 
distribution is not affected by these explosions. Only the hot 
gas is slightly removed from the central region of the galaxy 
by SNe. AGN feedback removes dense gas from the galaxy 
and pushes hot gas to very large distances. Furthermore the 
amount of cold filamentary material in the core of the halo 
(at 0.1 r v i r ) is decreased by one order of magnitude, suggest- 
ing that filaments are also pushed away to larger distances. 
This is confirmed by Fig. 10 showing an outflow of cold fila- 
ments at 7Vi r around z = 6. The removal of cold streams from 
the centre of the halo can potentially reduce the Lyman-a 
emission of simulated Lyman-a blobs, which usually tend 
to be brighter than observed (see Rosdahl & Blaizot 2012, 
where no AGN feedback is accounted for). However a full 
treatment of the radiative transfer from stars and the cen- 
tral AGN is required to track this emission. 

As visual inspection of the gas density in Fig. 7 al- 
ready suggests, the cold filamentary structures exhibit mor- 
phological differences between the nofeedback / no AGN cases 
and the AGN case. We quantify these differences using the 
skeleton (Sousbie et al. 2009), a geometric 3D ridge extrac- 
tor well suited to identify filaments. We build the skeleton 
of the gas density within the virial radius of the halo for 
the no AGN and the AGN cases at different redshifts, and 
compare the properties of the skeleton such as its average 
length and curvature. The gas density is sampled over a uni- 
form grid of 256 3 pixels inside a box of 1.1 comoving Mpc 
(4.5 r v ir = 160 physical kpc at z = 6) and is smoothed with 
a Gaussian filter of 4a. Once the skeleton is computed, it 
is thresholded both within r v i r and where gas number den- 
sity is larger than nn > 50 Hem -3 . The two skeletons are 
globally very similar but depart significantly in shape on 
small scales (Fig. 12), with the AGN skeleton being both 
shorter and curlier with decreasing redshift relative to the 
noAGN skeleton. More specifically, Fig. 12 shows the ratio 
between the AGN and the no AGN case of the total length 
of the thresholded skeletons as well as the ratio of their 
mean curvature. As expected, compared to the noAGN case, 
the total length of the skeleton in the AGN case decreases 
by about 50 per cent, while its curvature (as defined in 
Pogosyan et al. 2009) increases by about 20 per cent. These 
are clear signatures that AGN feedback is both destroying 
the connectivity of the filaments and bending them. 



6 CONCLUSION AND DISCUSSION 

Using high resolution cosmological hydrodynamical simula- 
tions to follow the formation of the progenitor of a BCG 
at z = 6, we have studied the impact of early BH growth 
and its accompanying AGN activity on the formation of its 
host proto-galaxy and on its environment. Our BH and AGN 
model is calibrated to match the BH-host relationships at 
z=0. The main results of this investigation are the following: 

• SN feedback is extremely inefficient at regulating star 
formation in progenitors of massive galaxies and has almost 
no impact on the surrounding gas of the galaxy. 
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Figure 10. Mass fluxes measured at r v i r as a function of redshift for the nofeedback case (black) noAGN case (blue), and for the 
AGN case (red). Net fluxes are solid lines, outflow components are dotted lines, and inflow components are dashed lines. Top left panel 
is the total flux of gas, top right panel are the cold flows and star forming gas, bottom left panel is the cold diffuse only, and bottom right 
is the hot phase only. Bottom panels have a different scaling of the y-axis than the top panels. Global accretion is significantly reduced 
due to AGN feedback, CD accretion is quenched, HD accretion is strongly reduced with a strong outflowing component, and the F+ISM 
component is slightly reduced with some cold gas being transported outside of the halo. 
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Figure 11. Integrated mass profiles at z = 6 of the cold filaments 
(solid), star forming gas (dotted), and hot gas (dashed) for the 
nofeedback (black), no AGN (blue), and AGN case (red). The ver- 
tical grey line corresponds to 0.1 r v i r . Profiles in the no AGN case 
are very similar to the nofeedback case: only a small quantity of 
hot gas is removed by SNe. AGN feedback however expels large 
amounts of gas from the centre of the halo. 



• AGN feedback through a quasar mode quenches the SF 
in the bulge and reduces the amount of gas accreted onto 
the central galaxy and available for star formation. 

• With rb = 50 pc and a stellar mass of Mb = 6.2 x 
10 9 Mq , the stellar bulge in our simulation with AGN feed- 
back is still a factor ten smaller than the most compact 
observed stellar bulges at 1 < z < 3 of comparable stellar 
mass. This suggests that the strong evolution of bulge sizes 
between z = and z = 2 continues to higher redshift. 

• The distribution of mass within the core of the halo is 
strongly affected by the central AGN: gas, stars and DM ex- 
hibit shallower density profiles in their centre. The adiabatic 
contraction of DM in the core of the halo is reduced by 35 
per cent due to the presence of the central AGN compared 
to the simulation without any feedback. 

• The cold baryon content is regulated by two stages of 
BH growth: (1) the rapid early growth of the BH at the Ed- 
dington limit driven by direct cold filamentary infall onto the 
bulge, followed by (2) a sub-Eddington phase where bursts 
of AGN activity are triggered by inflows and rapid clump 
migration towards the BH, brought on by merger induced 
torques. 
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Figure 12. Top: Image of the skeleton for the noAGN case (from 
black to green) and the AGN case (from black to orange). Bottom: 
The ratio of the AGN to the noAGN length (bottom bundle, red 
points) and curvature (top bundle, orange points) of the skele- 
ton within the halo as a function of redshift (shaded areas are 1 
and 2 a deviations). AGN feedback changes the structure of the 
filamentary network by bending and shortening it. 



• The self-regulation of the BH growth due to AGN feed- 
back limits the BH mass to a value a factor four below the 
z — MBH-fb relation and a factor 12 above the z = 
MeH-Afb relation. Our simulations therefore suggest that 
these relations evolve with redshift. 

• AGN feedback drives large-scale outflows that reduce 
the baryon fraction in the halo by 30 per cent. The gas that 
is expelled from the halo is essentially a hot component with 
temperature well above virial, and with an outflow velocity 
of 500 kms -1 at r vlI that is larger than the escape velocity 
of the halo. 

• Cold filaments are morphologically disturbed by AGN 
with a small fraction of them blown out of the halo, and a 
deficit of them in the core of the halo close to the galaxy at 
0.1r vir . 

For "canonical" feedback parameters (see section 2), 
we demonstrated the key role played by AGN feedback in 
regulating both the star formation, and the amount of gas 
available in the halo for the progenitor of a massive cluster 
of galaxies. We also tested two variations of our canonical 
model for AGN feedback with the following results. A model 



where AGN feedback energy is only deposited in a radio 
(jet) mode produces bipolar outflows similar to those pro- 
duced by our canonical model in which energy is released 
in isotropic thermal bursts; thus, the impact on the sur- 
rounding gas is as efficient as in the canonical case (see Ap- 
pendix A). Reducing the AGN feedback efficiency of our 
canonical model by one order of magnitude also does not 
change the AGN's impact on the star formation in the cen- 
tral galaxy: the BH grows to larger masses so as to release 
the same amount of energy to unbind the bulge gas (see 
Appendix B). These very simple tests lend support to our 
results on the impact of AGN feedback, showing them to be 
robust to variations of its modeling. In the future, alternative 
variations could also be explored: for instance, gravitational 
recoil of BHs which commonly takes place in massive ha- 
los (Volonteri et al. 2010); a mode of quasar heating where 
the radiative transfer of photons and their radiation pressure 
is followed self-consistently (Kim et al. 2011). 

Our numerical experiments have possible implications 
for models of pre-heating of the gas at high redshift. Stan- 
dard pre-heating models suppose that gas elements are 
more or less uniformly heated (either by imposing a tem- 
perature or an entropy floor for the gas) to some value 
large enough to prevent the formation of excessively lu- 
minous clusters and to break the self-similarity of clus- 
ters (Bialek et al. 2001; Babul et al. 2002; Borgani et al. 
2002; Muanwong et al. 2002; Voit et al. 2003). However, 
these models fail to predict some of the basic properties 
of clusters, producing large isentropic cores (Borgani et al. 
2005; Younger & Bryan 2007) and shallow pressure pro- 
files (Kay et al. 2012). Also the uniform pre-heating leaves 
the IGM with the wrong covering fraction for the Lyman- 
q forest unless only high gas density regions are pre- 
heated (Shang et al. 2007; Borgani & Viel 2009). Our simu- 
lations demonstrate that AGN feedback is a potential strong 
source of pre-heating, but instead of uniformly changing the 
specific energy of the gas, it changes the temperature, and 
thus the entropy of the hot diffuse gas component, leaving 
the thermodynamical (but not dynamical) properties of the 
cold streams unchanged. Thus, one can get a more realistic 
model of pre-heating by only changing the temperature of 
the gas that is already hot (above 10 5 K). 

There are several potential worries concerning our re- 
sults and the probable consequences on galaxy formation 
at lower redshift. The SFR measured in the central galaxy 
at 2 = 6 is only 20Moyr~ 1 . A naive back-of-the enve- 
lope calculation predicts that, with this amount of SFR, 
the final BCG mass would be of the order of 10 M© , 
one order of magnitude below the expected mass for BCGs 
(10 12 M©, Moster et al. 2010). However, the central galaxy 
can build a significant fraction of its mass through merger 
events. With the SFRs of our AGN feedback case, more 
than 90 per cent of the stellar mass would need to be added 
by ex-situ star formation (as opposed to in-situ). Hydro 
simulations for the most massive structures show that this 
is plausible (Oser et al. 2010). Our result is in agreement 
with observations of z — 6 galaxies of comparable stel- 
lar mass (~ 10 10 Mq) that show a SFR of the same level 
20M Q yr- 1 (Stark et al. 2009), even though taking into ac- 
count correct dust extinction would increase the SFR by a 
factor 2 (Bouwens et al. 2012). Also, the strong outflows and 
self-regulation of the cold baryons produced by our AGN 
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model provides a suitable solution to suppress the SFR be- 
tween 4 < z < 7 in order to reproduce the plateau of the 
specific SFR above z > 2 (Weinmann et al. 2011), and to 
decrease the specific star formation rate that is required for 
the most massive galaxies (Kimm et al. 2012). 

Another potential tension is the low baryon fraction ob- 
tained in the AGN feedback simulation, i.e. /b = 0.11 at z = 
6. Observed clusters in this mass range M vir = 2 x 10 15 M 
at z = exhibit a universal baryon fraction (Gonzalez et al. 
2007; Giodini et al. 2009). However, this discrepancy is an 
issue only if it is assumed that the gas cannot re-collapse. 
Simulations of strong pre-heating of structures show that 
the gas fraction is reduced early on, and is followed by 
re-accretion of the hot expelled gas that raises up the 
baryon fraction at lower redshift (Faucher-Giguere et al. 
2011; Peirani et al. 2012). The capacity of AGN feedback 
to expel gas from halos is a key problem for precision cos- 
mology and future surveys, as large-scale simulations have 
shown that the presence of baryons and their associated 
physics modify by a few percent the measurements of the 
matter power spectrum on Mpc scales where possible sig- 
nature of modified gravity are sought (Guillet et al. 2010; 
van Daalen et al. 2011). 

In order to put stronger observational constraints on 
our current model of AGN feedback and confirm or inval- 
idate the plausibility of our theoretical predictions, simu- 
lations of more massive halos at z = 6 have to be car- 
ried out (in the spirit of Di Matteo et al. 2012). Another 
direction should also be explored, by simulating the for- 
mation of intermediate-redshift 2 = 3 — 1 galaxies where 
the direct imaging of galaxies is able to unveil the dumpi- 
ness of these gas-rich objects (Elmegreen et al. 2007), and 
put more severe constraints on the AGN-galaxy connec- 
tion (Bournaud et al. 2011). 



ACKNOWLEDGMENTS 

The simulations presented here were run on the DiRAC fa- 
cility jointly funded by STFC, the Large Facilities Capital 
Fund of BIS and the University of Oxford. This research is 
part of the Horizon-UK project. CP thanks T. Sousbie and 
C. Gay for skeleton related codes. YD and JS acknowledge 
support by the ERC advanced grant (Dark Matters). 



References 

Babul A., Balogh M. L., Lewis G. F., Poole G. B., 2002, 

MNRAS, 330, 329 
Begelman M. C, Volonteri M., Rees M. J., 2006, MNRAS, 

370, 289 

Bell E. F., Mcintosh D. H., Katz N., Weinberg M. D., 2003, 
ApJ Sup., 149, 289 

Bialek J. J., Evrard A. E., Mohr J. J., 2001, ApJ, 555, 597 

Binney J., 1977, ApJ, 215, 483 

Birnboim Y., Dekel A., 2003, MNRAS, 345, 349 

Blumenthal G. R., Faber S. M., Flores R., Primack J. R., 
1986, ApJ, 301, 27 

Boehringer H., Voges W., Fabian A. C, Edge A. C, Neu- 
mann D. M., 1993, MNRAS, 264, L25 

Bondi H., 1952, MNRAS, 112, 195 



Booth C. M., Schaye J., 2009, MNRAS, 398, 53 
Booth C. M., Schaye J., 2011, MNRAS, 413, 1158 
Borgani S., Finoguenov A., Kay S. T., Ponman T. J., 

Springel V., Tozzi P., Voit G. M., 2005, MNRAS, 361, 

233 

Borgani S., Governato F., Wadsley J., Menci N., Tozzi P., 
Quinn T., Stadel J., Lake G., 2002, MNRAS, 336, 409 

Borgani S., Viel M., 2009, MNRAS, 392, L26 

Bournaud F., Dekel A., Teyssier R., Cacciato M., Daddi 
E., Juneau S., Shankar F., 2011, ApJ Let, 741, L33 

Bournaud F. et al., 2011, ArXiv e-prints 

Bouwens R. J. et al., 2012, ApJ, 754, 83 

Bower R. G., Benson A. J., Malbon R., Helly J. C, Frenk 

C. S., Baugh C. M., Cole S., Lacey C. G., 2006, MNRAS, 
370, 645 

Brighenti F., Mathews W. G., 2006, ApJ, 643, 120 

Brooks A. M., Governato F., Quinn T., Brook C. B., Wad- 
sley J., 2009, ApJ, 694, 396 

Bruce V. A. et al., 2012, ArXiv e-prints 

Briiggen M., Ruszkowski M., Hallman E., 2005, ApJ, 630, 
740 

Cano-Di'az M., Maiolino R., Marconi A., Netzer H., Shem- 

mer O., Cresci G., 2012, A&A, 537, L8 
Cattaneo A., Dekel A., Devriendt J., Guiderdoni B., 

Blaizot J., 2006, MNRAS, 370, 1651 
Chapon D., Mayer L., Teyssier R., 2011, ArXiv e-prints 
Chartas G., Brandt W. N., Gallagher S. C, 2003, ApJ, 595, 

85 

Cimatti A. et al., 2008, A&A, 482, 21 
Croton D. J. et al., 2006, MNRAS, 365, 11 
Daddi E. et al., 2005, ApJ, 626, 680 

Decarli R., Falomo R., Treves A., Labita M., Kotilainen 

J. K., Scarpa R., 2010, MNRAS, 402, 2453 
DeGraf C, Di Matteo T., Khandai N., Croft R., 2012, ApJ 

Let, 755, L8 
Dekel A. et al., 2009, Nature, 457, 451 
Di Matteo T., Colberg J., Springel V., Hernquist L., Sijacki 

D. , 2008, ApJ, 676, 33 

Di Matteo T., Khandai N., DeGraf C, Feng Y., Croft 
R. A. C, Lopez J., Springel V., 2012, ApJ Let, 745, L29 

Di Matteo T., Springel V., Hernquist L., 2005, Nature, 433, 
604 

Dubois Y., Devriendt J., Slyz A., Teyssier R., 2010, MN- 
RAS, 409, 985 

Dubois Y., Devriendt J., Slyz A., Teyssier R., 2012, MN- 
RAS, 420, 2662 

Dubois Y., Devriendt J., Teyssier R., Slyz A., 2011, MN- 
RAS, 417, 1853 

Dubois Y., Pichon C, Haehnelt M., Kimm T., Slyz A., 
Devriendt J., Pogosyan D., 2012, MNRAS, 423, 3616 

Dubois Y., Teyssier R., 2008, A&A, 477, 79 

Elmegreen B. C, Bournaud F., Elmegreen D. M., 2008, 
ApJ, 688, 67 

Elmegreen D. M., Elmegreen B. G. Ravindranath S., Coe 

D. A., 2007, ApJ, 658, 763 
Fabjan D., Borgani S., Tornatore L., Saro A., Murante G., 

Dolag K., 2010, MNRAS, 401, 1670 
Fan X. et al., 2001, A J, 122, 2833 
Fan X. et al, 2006, AJ, 132, 117 

Faucher-Giguere C.-A., Keres D., Ma C.-P, 2011, MNRAS, 
417, 2982 



AGN and high-redshift galaxy formation 15 



Forman W., Kellogg E., Gursky H., Tananbaum H., Giac- 

coni R., 1972, ApJ, 178, 309 
Gaibler V., Khochfar S., Krause M., 2011, MNRAS, 411, 

155 

Gaspari M., Melioli C, Brighenti F., D'Ercole A., 2011, 

MNRAS, 411, 349 
Giodini S. et al., 2009, ApJ, 703, 982 

Gnedin O. Y., Kravtsov A. V., Klypin A. A., Nagai D., 

2004, ApJ, 616, 16 
Goerdt T., Dekel A., Sternberg A., Gnat O., Ceverino D., 

2012, MNRAS, 424, 2292 
Gonzalez A. H., Zaritsky D., Zabludoff A. I., 2007, ApJ, 

666, 147 

Guillet T., Teyssier R., Colombi S., 2010, MNRAS, 405, 
525 

Gursky H., Kellogg E., Murray S., Leong C, Tananbaum 

H., Giacconi R., 1971, ApJ Let, 167, L81 
Haardt F., Madau P., 1996, ApJ, 461, 20 
Raring N., Rix H.-W., 2004, ApJ Let, 604, L89 
Hopkins P. F., Hernquist L., Cox T. J., Di Matteo T., 

Robertson B., Springel V., 2006, ApJ Sup., 163, 1 
Jiang L. et al., 2009, AJ, 138, 305 

Kay S. T., Peel M. W., Short C J., Thomas P. A., Young 
O. E., Battye R. A., Liddle A. R., Pearce F. R., 2012, 
MNRAS, 422, 1999 

Keres D., Katz N., Fardal M., Dave R., Weinberg D. H., 
2009, MNRAS, 395, 160 

Keres D., Katz N., Weinberg D. H., Dave R., 2005, MN- 
RAS, 363, 2 

Khandai N., Feng Y., DeGraf C, Di Matteo T., Croft 

R. A. C, 2012, MNRAS, p. 3027 
Khochfar S., Silk J., 2006, ApJ Let, 648, L21 
Kim J.-h., Wise J. H., Alvarez M. A., Abel T., 2011, ApJ, 

738, 54 

Kimm T. et al., 2012, MNRAS, 425, L96 
Kimm T., Slyz A., Devriendt J., Pichon C, 2011, MNRAS, 
413, L51 

King A., 2003, ApJ Let., 596, L27 

Komatsu E. et al., 2011, ApJ Sup., 192, 18 

Krumholz M. R., Tan J. C, 2007, ApJ, 654, 304 

Leccardi A., Molendi S., 2008, ,40,4, 487, 461 

Li Y. et al, 2007, ApJ, 665, 187 

Magorrian J. et al., 1998, AJ, 115, 2285 

Maiolino R. et al., 2005, A&A, 440, L51 

Maiolino R. et al., 2012, MNRAS, 425, L66 

Markevitch M., Forman W. R., Sarazin C. L., Vikhlinin A., 
1998, ApJ, 503, 77 

Martizzi D., Teyssier R., Moore B., Wentz T., 2012, MN- 
RAS, 422, 3081 

McCarthy I. G., Schaye J., Bower R. G., Ponman T. J., 
Booth C. M., Dalla Vecchia C, Springel V., 2011, MN- 
RAS, 412, 1965 

McCarthy I. G. et al., 2010, MNRAS, 406, 822 

Meier D. L., Sadun A. C, Lind K. R., 1991, ApJ, 379, 141 

Merloni A. et al., 2010, ApJ, 708, 137 

Mortlock D. J. et al., 2011, Nature, 474, 616 

Moster B. P., Somerville R. S., Maulbetsch C, van den 
Bosch F. C, Maccio A. V., Naab T., Oser L., 2010, ApJ, 
710, 903 

Muanwong O., Thomas P. A., Kay S. T., Pearce F. R., 

2002, MNRAS, 336, 527 
Nagai D., Kravtsov A. V., Vikhlinin A., 2007, ApJ, 668, 1 



Ocvirk P., Pichon C, Teyssier R., 2008, MNRAS, 390, 1326 
Omma H., Binney J., Bryan G., Slyz A., 2004, MNRAS, 
348, 1105 

Oser L., Ostriker J. P., Naab T., Johansson P. H., Burkert 

A., 2010, ApJ, 725, 2312 
Ostriker E. C, 1999, ApJ, 513, 252 

Panter B., Jimenez R., Heavens A. F., Chariot S., 2007, 

MNRAS, 378, 1550 
Peirani S., Jung I., Silk J., Pichon C, 2012, ArXiv e-prints 
Pichon C, Pogosyan D., Kimm T., Slyz A., Devriendt J., 

Dubois Y., 2011, MNRAS, 418, 2493 
Pogosyan D., Pichon C, Gay C, Prunet S., Cardoso J. F., 

Sousbie T., Colombi S., 2009, MNRAS, 396, 635 
Pontzen A., Governato F., 2012, MNRAS, 421, 3464 
Powell L. C, Slyz A., Devriendt J., 2011, MNRAS, 414, 

3671 

Pratt G. W., Bohringer H., Croston J. H., Arnaud M., 
Borgani S., Finoguenov A., Temple R. F., 2007, A&A, 
461, 71 

Read J. I., Gilmore G., 2005, MNRAS, 356, 107 
Rosdahl J., Blaizot J., 2012, MNRAS, 423, 344 
Rosen A., Bregman J. N., 1995, ApJ, 440, 634 
Ruszkowski M., Briiggen M., Begelman M. C, 2004, ApJ, 
611, 158 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Shang C, Crotts A., Haiman Z., 2007, ApJ, 671, 136 
Sijacki D., Springel V., Di Matteo T., Hernquist L., 2007, 

MNRAS, 380, 877 
Sijacki D., Springel V., Haehnelt M. G., 2009, MNRAS, 

400, 100 

Silk J., Rees M. J., 1998, A&A, 331, LI 

Somerville R. S., Hopkins P. F., Cox T. J., Robertson B. E., 

Hernquist L., 2008, MNRAS, 391, 481 
Sousbie T., Colombi S., Pichon C, 2009, MNRAS, 393, 457 
Springel V., Di Matteo T., Hernquist L., 2005, MNRAS, 

361, 776 

Stark D. P., Ellis R. S., Bunker A., Bundy K., Targett T., 

Benson A., Lacy M., 2009, ApJ, 697, 1493 
Sternberg A., Pizzolato F., Soker N., 2007, ApJ Let, 656, 

L5 

Sutherland R. S., Bicknell G. V., 2007, ApJ Sup., 173, 37 
Sutherland R. S., Dopita M. A., 1993, ApJ Sup., 88, 253 
Teyssier R., 2002, A&A, 385, 337 

Teyssier R., Chapon D., Bournaud F., 2010, ApJ Let, 720, 
L149 

Teyssier R., Moore B., Martizzi D., Dubois Y., Mayer L., 

2011, MNRAS, 414, 195 
Tornatore L., Borgani S., Springel V., Matteucci F., Menci 

N., Murante G., 2003, MNRAS, 342, 1025 
Treister E., Urry C. M., Schawinski K., Cardamone C. N., 

Sanders D. B., 2010, ApJ Let, 722, L238 
Tremaine S. et al., 2002, ApJ, 574, 740 
Valdarnini R., 2003, MNRAS, 339, 1117 
van Daalen M. P., Schaye J., Booth C. M., Dalla Vecchia 

C, 2011, MNRAS, 415, 3649 
Vikhlinin A., Markevitch M., Murray S. S., Jones C, For- 
man W., Van Speybroeck L., 2005, ApJ, 628, 655 
Voit G. M., Balogh M. L., Bower R. G., Lacey C. G., Bryan 

G. L., 2003, ApJ, 593, 272 
Volonteri M., Giiltekin K., Dotti M., 2010, MNRAS, 404, 

2143 

Volonteri M., Rees M. J., 2006, ApJ, 650, 669 



16 Y. Dubois et al. 




Figure Al. Slices of the gas density at z = 7.5 for the AGN run 
where the energy is released through thermal bursts (top panel), 
or through collimated jets (bottom panel). Arrows indicate the ve- 
locity vectors in the plane of the image, and their lengths are pro- 
portional to the velocity amplitudes. The image sizes are 15 kpc. 
In both cases, an expanding wind develops from the central disc. 
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APPENDIX A: ISOTROPIC WIND VERSUS 
COLLIMATED JET INPUTS 

Our standard implementation of AGN feedback deposits en- 
ergy in isotropic thermal bursts when BHs accrete at high 



rates (x 0.01) and collimated jets at low accretion rates 
(x < 0.01). Both the thermal inputs and the jets deposit 
their energy (mass and momentum for the jet) in a small re- 
gion that is approximated by a sphere of radius Ax = 15 pc 
for the quasar mode, and a cylinder of radius Ax and height 
2 Ax for the radio mode. The size of the region of the energy 
input is typically two orders of magnitude smaller than the 
characteristic radius of the galactic disc (2 kpc at z = 6). In 
the radio mode, the momentum injected in each cell within 
the cylinder is parallel to the axis defined by the angular mo- 
mentum vector of the gas surrounding the BH (maximum 
collimation with no opening angle). As shown in Fig. 2, the 
BH accretes gas above x ^ 0.01 all the way down to redshift 
z — 6, thus energy is released through thermal pulses that 
should propagate isotropically in a homogeneous medium. 
However, the BH is surrounded by a thin disc of dense star- 
forming gas. As a result the outflow produced by the AGN 
in its quasar mode is not perfectly isotropic but exhibits 
a bipolar outflow with a large opening angle (Fig. Al). It 
is possible that in reality winds are much more collimated, 
and that our predictions for the impact of AGN feedback 
on the galaxy and its surrounding gas are overestimated be- 
cause the wind has a cross-section that is too large. To test 
this, we run a simulation where AGN feedback is only al- 
lowed to release its energy in a collimated jet mode (the 
so-called radio mode). The result is illustrated in the bot- 
tom panel of Fig. Al. The large-scale outflow produced by 
this jet mode is very similar to the isotropic thermal wind: a 
bipolar expanding wind develops around the central galaxy. 
The reason for the quick expansion of the radio AGN is that 
the jet is light. As a consequence the pressure inside the jet 
pushes the cavity surrounding the jet and the wind expands 
as much as it can in all directions (Meier et al. 1991) The 
multiphase structure of the interstellar medium also help to 
destabilize the jet propagation (Sutherland & Bicknell 2007; 
Gaibler et al. 2011). 



APPENDIX B: REDUCING AGN FEEDBACK 
EFFICIENCY 

We also test how the growth of the BH and the impact of 
AGN feedback when the efficiency £f of AGN energy cou- 
pling to the gas is reduced by a factor of 15. In this case the 
central supermassive BH reaches a larger mass Mbh.Iow = 
5 x 10 s M at z = 6.7 instead of M BH ,Wgh = 6 x 10 7 M 
for our canonical model at the same redshift. Thus, the en- 
ergy released by the two BHs are very similar (the energy 
released by the canonical model is 1.8 times the energy of 
model with reduced efficiency) because the energy required 
to unbind the dense cold gas surrounding the BH does not 
strongly vary with time. BHs accrete close to, or at their Ed- 
dington rate before blowing out the gas and self-regulating 
their growth. The characteristic growth time of a BH ac- 
creting at the Eddington limit is t^dd — 45 Myr for a ra- 
diative efficiency of e r = 0.1. Thus, two BHs accreting gas 
at the Eddington limit with AGN feedback efficiencies e t 
that are one order of magnitude different will reach the self- 
regulated state with a time delay of the order of iEdd (see 
Fig. Bl). As the AGN feedback luminosities are compara- 
ble in both runs, the bolometric luminosity is naturally one 
order of magnitude larger in the simulation with reduced effi- 
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Figure Bl. AGN feedback luminosity of the most massive BH 
for the canonical AGN model (black), the reduced efficiency re- 
duced by 15 (grey), and their bolometric luminosities (red and 
orange respectively). The feedback luminosities are similar for 
the two AGN feedback models, but their bolometric luminosities 
arc different by one order of magnitude. 
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Figure B2. SFR history using the stars contained within 0.1 r w \ T 
for the nofeedback case (black), the no AGN case (blue), and for 
the AGN case with the canonical AGN feedback model (red), 
and the model with reduced efficiency (orange). The two AGN 
feedback models show similar SFRs with values reduced by a 
factor of 10 compared to nofeedback /no AGN cases. 



APPENDIX C: RESOLUTION TESTS 

We added some other simulations where the minimum al- 
lowed cell size is degraded down to Ax = 135 pc (above 
the nominal resolution used in the paper of Ax = 15 pc) 
to test the influence of resolution on the impact of AGN 
feedback. One of these simulations is the equivalent of the 
nofeedback case (without SN and no AGN feedback), and 
another one is the equivalent of the AGN case (with SN 
and AGN feedback) described in section 2. The change in 
resolution is accompanied by a modification in the threshold 
for star formation, which is now po = 5 Hem -3 . The con- 
sequence is that the star forming gas is more supported by 
pressure due to a lower value of po in the polytropic equa- 
tion of state and produces a less clumpy medium with the 
star formation taking place more smoothly. We added one 
last run by degrading both the DM resolution and spatial 
resolution to M rcs = 7 x 10 6 h' 1 M and Ax = 135 pc with 
the same physics as our canonical AGN case, and for which 
the threshold for star formation is po = 5Hcm~ 3 . 

Fig. CI shows the SFRs measured within O.lrvir for 
these low resolution simulations and compared to the high 
resolution simulations. The SFRs with the same physics 
but different resolutions have very similar values between 
6 < z < 7, which suggests that our results are well converged 
in this redshift range: the starvation of the gas on the cen- 
tral galaxy by AGN quenching has a comparable impact on 
the SFR at late times. Above z = 7, the SFRs are different 
from on resolution to another. We interpret this as a result 
of both the force being sampled on coarser meshes, delaying 
the collapse of structures, and of a less clumpy and less dense 
ISM at low resolution, and as a consequence, the time-scales 
associated to star formation are longer (see Teyssier et al. 
2010). The measured baryon fraction at z — 6 and at T*vir IS 
reduced by 30 per cent below the universal value for low res- 
olution and high resolution AGN simulations. We also find 
that, at low resolution, the fraction of cold flows is reduced 
by a factor 2 at 0.1r v ir due to the presence of AGN feedback 
compared to the lo-resolution nofeedback case (at high res- 
olution this reduction factor is 3, see Fig. 11). We can infer 
that the results found in this paper are robust to a dramatic 
change (one order of magnitude) in DM mass resolution and 
in minimum cell size. 



ciency, and reaches values up to 10 L© compatible with the 
bolometric luminosities observed for the brightest quasars at 
z ~ 6 (10 12 - 10 14 L Q , Willott et al. 2010). As the energy 
released by the central BHs is comparable, the impact on 
the galaxy star formation is also similar (see Fig. B2). The 
two models of AGN feedback significantly reduce the SFR 
in the central object with levels of star formation that are 
comparable. Hence the main results of this paper do not de- 
pend on the details of the feedback implementation and is 
only weakly dependent on the feedback efficiency 6f . 
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Figure CI. SFR history using the stars contained within 0.1 r v i r 
for the nofeedback case with Ax = 15 pc (black dashed), the 
AGN case with Ax = 15 pc (red dashed), nofeedback case with 
Ax = 135 pc (black solid), the AGN case with Ax = 135 pc 
(red solid), and the AGN case with Ax = 135 pc and M TeB = 
7 X 10 6 h~ x Mq (green). The SFRs with different resolution using 
the same physics have similar values at z = 6, but they show a 
significant difference at higher rcdshift z > 7. 



